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Co-doping of Ce-doped LaBra with Ba, Ca, or Sr improves the energy resolntion that can be 
achieved by radiation detectors based on these materials. Here, we present a mechanism that 
rationalizes of this enhancement that on the basis of hrst principles electronic structnre calculations 
and point defect thermodynamics. It is shown that incorporation of Sr creates neutral Vbi — Srca 
complexes that can temporarily trap electrons. As a result, Auger quenching of free carriers is 
reduced, allowing for a more linear, albeit slower, scintillation light yield response. Experimental 
Stokes shifts can be related to different Ceua — Srca — Vbv triple complex configurations. Co-doping 
with other alkaline as well as alkaline earth metals is considered as well. Alkaline elements are found 
to have extremely small solubilities on the order of 0.1 ppm and below at 1000 K. Among the alkaline 
earth metals the lighter dopant atoms prefer interstitial-like positions and create strong scattering 
centers, which has a detrimental impact on carrier mobilities. Only the heavier alkaline earth 
elements (Ca, Sr, Ba) combine matching ionic radii with sufficiently high solubilities. This provides 
a rationale for the experimental Ending that improved scintillator performance is exclusively achieved 
using Sr, Ca, or Ba. The present mechanism demonstrates that co-doping of wide gap materials 
can provide an efficient means for managing charge carrier populations under out-of-equilibrium 
conditions. In the present case dopants are introduced that manipulate not only the concentrations 
but the electronic properties of intrinsic defects without introducing additional gap levels. This 
leads to the availability of shallow electron traps that can temporarily localize charge carriers, 
effectively deactivating carrier-carrier recombination channels. The principles of this mechanism 
are therefore not specific to the material considered here but can be adapted for controlling charge 
carrier populations and recombination in other wide gap materials. 

PACS numbers: 61.72.-y 61.72.Bb 72.20.Jv 78.70.Ps 


I. INTRODUCTION 

Many applications in nuclear and radiological surveil¬ 
lance, high-energy physics, and medical imaging rely on 
scintillator materials,which enable the energy resolved 
detection of high energy radiation.^ The energy resolu¬ 
tion that can be accomplished increases with luminosity, 
which is usually related to the efficiency of the process by 
which the energy of incoming radiation quanta (typically 
in the keV to MeV range) is converted to lower energy 
photons (on the order of a few eV). The achievable res¬ 
olution is, however, further limited by the non-linearity 
of the scintillation response to the energy of the incident 
radiation,^ which arises from the competition between 
non-radiative quenching, defect carrier trapping, as well 
as activator capture and subsequent emission.®^® 

Recent work showed that the energy resolution of Ce- 
doped LaBra can be significantly improved by co-doping. 
The concept was first realized experimentally by Yang et 
al. for samples co-doped with Sr.^° Later Alekhin et al. 
achieved an improvement of energy resolution from 2.7 
to 2.0% at 662 keV using Sr and Ca.^^ A more compre¬ 
hensive investigation including both the alkaline as well 
as alkaline earth series revealed that better performance 
is only achieved when using the heavier elements of the 
latter series (Sr, Ca, Ba).^^ To explain these observa¬ 
tions it has been suggested that doping with Sr, Ca, or 
Ba causes^^ (*) s- reduction of the non-radiative recombi¬ 


nation rate, {ii) an increase of the so-called escape rate 
of the carriers from the quenching region, or {Hi) an in¬ 
crease in the trapping rate of Ce^^. The experimental 
investigations also revealed three distinct optical signa¬ 
tures associated with Ce, which have been interpreted as 
evidence for the presence of three different Ce environ¬ 
ments in the co-doped material. In contrast, only one 
such feature can be identified in LaBr 3 :Ce. Later, it was 
argued by the present authors that Sr-doping causes the 
creation of shallow electron trap complexes, which leads 
to reduced Auger quenching.^"* 

The present paper describes the argumentation in de¬ 
tail and presents a careful analysis of the thermody¬ 
namic properties and electronic structures of the most 
important intrinsic and extrinsic defects - including their 
complexes - in Ce and Sr-doped LaBra including self- 
trapped polaronic configurations. To demonstrate that 
the present model is consistent with experimental obser¬ 
vations, a comprehensive set of calculated absorption and 
emission energies for the relevant Ce complexes is carried 
out. We obtain very good overall agreement, which en¬ 
ables us to correlate optical signatures with individual 
defect configurations. 

Finally, it is shown that the solubilities of Sr, Ca, and 
Ba near the synthesis temperature are several 100 ppm, 
whereas much lower values are obtained for the alkaline 
metals. Be and Mg yield large solubilities but also cause 
large lattice distortions and effectively act as interstitials. 
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These findings provide a rationale for why only Sr, Ca, 
and Ba have been experimentally found to improve scin¬ 
tillation performance. 

As part of the present work we also introduce a con¬ 
venient method for generating supercells with optimal 
shapes for defect calculations and finite-size scaling. The 
approach described in the appendix is applicable to ar¬ 
bitrary lattice types. 

The paper is organized as follows. Section II describes 
the computational methodology used in this work. The 
results regarding the thermodynamic properties of intrin¬ 
sic defects as well as Sr and Ce dopants are presented 
in Sect. III. Stokes shifts for single substitutional Ce 
and the triple complexes are reviewed and compared to 
experimental data and a dopant solubility analysis are 
presented. To conclude, we discuss our findings in light 
of scintillator performance in Sect. IV. The Appendices 
provide additional information concerning defect thermo¬ 
dynamics, solubilities, and finite size scaling of formation 
energies. 


II. METHODOLOGY 


A. General computational parameters 


Calculations were performed within density functional 
theory (DFT) using the projector augmented wave 
(PAW) method^^ as implemented in the Vienna ab- 
initio simulation package.^® Exchange-correlation (XC) 
effects were treated within the generalized gradient spin 
approximation.^^ DFT-l-17 type on-site potentials^® were 
included for both La-4/ {UeS = 10.3 eV) and Ce-4/ states 
{Ueff = 1.2 eV) in order to obtain the correct ordering of 
La-5(i and 4/ states and to reproduce experimental Ce- 
4/ ionization energies.The plane wave energy cutoff 
was set to 230 eV and Gaussian smearing with a width 
of 0.1 eV was used to determine occupation numbers. 

Lanthanum bromide adopts a hexagonal lattice struc¬ 
ture in space group 176 (POs/m) with La and Br ions 
occupying Wyckoff sites 2c and Qh, respectively. Using 
the aforementioned computational parameters one ob¬ 
tains lattice parameters of a = 8.140 A and c = 4.565 A to 
be compared with experimental values of a = 7.9648(5) A 
and c = 4.5119(5) A, respectively.^^ 

Spin-orbit interaction was not included self- 
consistently but rather added as a perturbation to 
the 4/ states for optical transitions according to 


AEs, 


-2Uf J=5/2 

V 2 C 4 / j = V 2 , 


where ^ 4 / = 0.1 eV as obtained from the 4/ splitting in 
a Ce-4/®5d^ configuration. 


B. Excited states 

Excited Ce-4/® states were obtained in a similar ap¬ 
proach as used by Canning and co-workers.First, a 
subspace of Ce-4/ states is identified by selecting states 
whose projection onto spherical harmonics with I = 3 
exceeds a suitably chosen threshold. This is made pos¬ 
sible by the localized nature of the rare-earth 4/ states. 
In this particular case, a Ce 4/^5(i® —)■ 4/®5fi^ excita¬ 
tion was emulated by completely deoccupying the 4/ 
manifold while keeping the number of electrons fixed. 
This results in a narrow, almost degenerate, group of 
unoccupied group of 14 4/ bands inside the fundamental 
electronic gap. The occupation numbers of Kohn-Sham 
states above this group is then occupied using Gaussian 
smearing, allowing for excited state structural relaxation. 


C. Defect calculations 

Defect formation energies were calculated using a well 
established thermodynamic formalism^"^ that is summa¬ 
rized for convenience in Appendix A. Lanthanum bro¬ 
mide has a small dielectric constant and some defects 
adopt large charge states (|g| < 4). As a result, im¬ 
age charge interactions are substantial and care must be 
taken to remove finite-size effects in order to obtain for¬ 
mation energies for the dilute limit.To this end a care¬ 
ful finite-size scaling study was carried out that is sum¬ 
marized in Appendix C. This Appendix also introduces 
a general and convenient method for constructing suit¬ 
able supercells with optimal shape given a certain sys¬ 
tem size, which allows one to obtain a dense sampling 
of different system sizes also in the case of low crystal 
symmetry. Hole polarons were studied using the polaron 
self-interaction correction method (pSIC).^® 

The bulk of the data presented in the following were 
obtained using 168-atom supercells. F-point sampling 
was found to be sufficient to converge defect formation 
energies to better than 0.05 eV. Configurations were re¬ 
laxed until ionic forces were less than lOmeV/A. Poten¬ 
tial alignment as well as periodic image charge correc¬ 
tions were taken into account to correct for finite size 
effects as detailed in Appendix C.^® 

It is well known that conventional XC functionals in¬ 
cluding the ones used in the present work commonly lead 
to an underestimation of band gaps. As discussed for ex¬ 
ample in Refs. 24, 25, 27-29, this error also affects defect 
formation energies and therefore also defect concentra¬ 
tions. Hybrid XC functionals, which combine conven¬ 
tional DFT functionals with Fock exchange, are often 
found to improve the band gap and are therefore ex¬ 
pected to also yield improved formation energies. For the 
present case, however, hybrid functionals are ill-suited 
since they cannot even qualitatively describe the posi¬ 
tion of the occupied Ce-4/ level in LaBra. This is re¬ 
lated to the distinct character of the electronic states in¬ 
volved. They are associated with widely different levels 
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of localization and thus the effective screening cannot be 
parametrized using a single (static) mixing parameter. 

In the present study we therefore resort to a simple 
correction scheme^®’^^^^® that shifts the formation ener¬ 
gies based on the offset between the “true” band edges 
and the ones obtained within the underlying computa¬ 
tional framework, i.e., in the present case DFT-|-[/ cal¬ 
culations. The offsets were determined using the Golfb 
method,which was previously found to yield a much 
improved description of the band structure of LaBrs com¬ 
pared to both DFT and DFT-|-{7.^^ Calculations for the 
primitive cell were carried out using a F-centered 3x3x6 
fc-point mesh and PAW data sets optimized for GW cal¬ 
culations, which require also unoccupied higher energy 
states to be well described. The dielectric tensor was 
computed for energies up to 200 eV above the CBM, 
equivalent to 1024 unoccupied bands. The offsets ob¬ 
tained in this way are AAvb = — 1.2eV for the valence 
band edge and AEcb = +0.5 eV for the conduction band 
edge. This increases the DFT+f7 band gap from 3.6 eV 
to 5.3 eV, which is in much better agreement with the ex¬ 
perimental value of 5.9 eV.^^ The correction scheme was 
only applied to defect charge states that did not include 
occupied localized states inside the band gap. As this 
distinction can be ambiguous, additional GoWo calcula¬ 
tions were carried out for 96-atom supercells of the most 
important defect configurations to verify the results from 
the correction scheme. In these calculations the Brillouin 
zone was sampled using the F-point only and the same 
GIF optimized PAW data sets as before were employed. 
The dielectric function was calculated up to 36 eV above 
the CBM. Based on this comparison the error in the for¬ 
mation energies is estimated to be 0.2 eV or less whereas 
for the transition levels the error is estimated to be below 
0.1 eV. 

Finally, defect concentrations were obtained using the 
calculated formation energies on the basis of a self- 
consistent solution of the charge neutrality condition, 
which has been described in detail in Refs. 24 and 29. 


III. RESULTS 


A. Intrinsic defects 


In this section the defect equilibria in nominally pure 
LaBra are investigated. To this end, one first requires 
knowledge of the formation energies of intrinsic defects, 
which are shown for La and Br-rich conditions (see Ap¬ 
pendix A) in Fig. l(a,b). It is apparent that under both 
La and Br-rich conditions the most important intrinsic 
donor and acceptor defects are Br and La vacancies, re¬ 
spectively, with interstitials and antisites playing minor 
roles. 

Given the respective formation energies it is straight¬ 
forward to compute the equilibrium transition level be¬ 


tween charge states qi and q 2 according to 

(?2 - <?1 


( 1 ) 


Experimentally, these can be detected for example by 
deep level transient spectroscopy. The levels calculated 
according to Eq. (1) are shown as red bars (left hand 
columns) in Fig. 3, which reveals that both types of va¬ 
cancies are associated with deep equilibrium transition 
levels. Note in particular that the +1/0 equilibrium tran¬ 
sition level of the Br vacancy is located 0.55 eV below the 
CBM. In GqIFo calculations, the quasi-particle energies 
associated with Vbv defect levels shift by less than 0.1 eV 
relative to the CBM, which provides strong support for 
this positioning of the defect level with respect to the 
conduction band edge. This is illustrated in Fig. 2(c,d), 
which compares quasi-particle energies from DFT and 
GqIFo calculations for both isolated and complexes Br 
vacancies. Note that for consistency, only transition lev¬ 
els obtained from DFT+I/ calculations including band 
gap corrections are reported from here on. 

Deep defects are typically associated with pronounced 
changes in the ionic positions between different charge 
states. This is indeed shown to be the case for Br va¬ 
cancies in Fig. 2(a), which shows the potential energy 
surfaces (PES) for both neutral and charged vacancies 
along the configuration coordinate connecting the respec¬ 
tive minima. A particular configuration is given in terms 
of the configuration coordinate ^ as 

-R = - Rtiln) /«FC, (2) 

where R^in denotes the minimum of the PES for charge 
state q and opc = \Rmin~Rtvin \ measures the structural 
difference between the two geometries. 

In the ideal structure, Br sites are surrounded by three 
cations (two out-of-plane and one in-plane with respect 
to the {0001} basal plane) at distances between 3.1 and 
3.2 A and eight Br ions (two in-plane, six out-of-plane) 
at distances of 3.6-3.7 A, see Fig. 4. In the charged va¬ 
cancy conhguration, the La neighbors of the vacant site 
move outward by 0.2-0.3 A while the Br neighbors move 
inward by up to 0.3 A. These relaxations are inverted for 
the neutral vacancy as La and Br neighbors are shifted 
inward and outward, respectively compared to the ideal 
structure. The large differences in ionic configuration are, 
as shown in Fig. 2(a), associated with substantial relax¬ 
ation energies of 0.25 and 0.38 eV on the neutral and 
charged PES, respectively. As illustrated in Fig. 4(c), 
the defect level is predominantly composed of 5d states 
localized at the three cations surrounding the vacant Br 
site, which as discussed below is of crucial importance 
for understanding the effect of Br vacancies on the op¬ 
tical signature of Ce. The defect level can act as an 
efficient electron trap, effectively removing carriers from 
the light-generation process during the instrumentation 
pulse shape-time. 

In terms of electronic trapping, also “vertical” transi¬ 
tions are important, which are indicated by the vertical 
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FIG. 1. Defect formation energies for (a,b) intrinsic defects including Vk-centers as well as (c,d) Sr. Formation energies for Sr 
related defects where computed assuming equilibrium with SrBr 2 (compare Sect. HIE). 
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FIG. 2. (a) Configuration coordinate diagram illustrating the relation between neutral and charged Br vacancies. For conve¬ 
nience the electron chemical potential was assumed to coincide with the CBM (^e = ccbm) such that the energy differences 
indicated by the vertical arrows correspond directly to the transition levels indicated in the middle panel, (b) Formation 
energy of Br vacancy configurations for La-rich conditions as a function of the electron chemical potential in the vicinity of the 
conduction band edge, compare Fig. 1(a). The solid and dotted lines are based on fully relax;ed configurations, corresponding 
to minima in panel (a). Dashed lines indicate formation energies representing vertical transitions in panel (a). The horizontal 
blue and orange arrows thus correspond to the vertical arrows in panel (a). The red horizontal arrow is equivalent to the energy 
difference between the minima in panel (a). (c,d) Quasi-particle spectra from DFT and GqWo calculations for both (c) the 
isolated Br vacancy and (d) a SrLa — Vbt complex. The position of the defect level relative to the CBM changes by less than 
0.1 eV when comparing the two types of calculations. Note that in both cases the defect level is unoccupied in equilibrium. 


arrows in Fig. 2(a). They can be calculated in a fashion 
similar to Eq. (1) but imposing the constraint of fixed 
ionic position. Specifically, the trapping level for charge 
state qi is obtained as 

etr( 9 i/ 92 ) =-, (o) 

92 - 9i 

where AEf{q 2 ; is the formation energy computed 

for charge state 92 at the ionic coordinates corresponding 
to the equilibrium positions in charge state 91 . All trap¬ 
ping levels presented in this paper were calculated with 
respect to the respectively more favorable charge state, 
e.g., 9 = -1-1 —)■ 0 for the Br vacancy and —3 —)■ —2, 


—2 —>■ — 1 ... for the La vacancy. 

Equilibrium transition and trapping levels can also be 
identified with crossing points in plots of the formation 
energy vs. electron chemical potential as illustrated in 
Eig. 2(b). In this figure, dashed lines correspond to for¬ 
mation energies computed in the geometry of the respec¬ 
tive other charge state, i.e. in the case of the neutral 
vacancy at versa. The equivalence of 

trapping and detrapping type transitions relative to the 
charge vacancy are shown by the blue and orange arrows 
in Figs. 2(a) and (b). Both equilibrium transition and 
trapping levels are compiled in Fig. 3 for the most rele¬ 
vant intrinsic and extrinsic (see below) defects. 
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FIG. 3. Equilibrium transition in red (left hand columns) and 
trapping levels in blue (right hand columns) obtained in the 
fashion indicated in Fig. 2. 



FIG. 4. Structure of lanthanum bromide and defect conhgu- 
rations therein, (a) Projection onto the (0001) basal plane. 
Large blue and small red spheres represent La and Br sites, re¬ 
spectively. The shaded triangles correspond to atomic layers 
containing both atom types (green) and Br ions only (lila), re¬ 
spectively. (b) Perspective view. A representative Br vacancy 
site is indicated by the cube and Sr as well as Ge substitu¬ 
tional sites are shown by large green and light gray spheres, 
respectively, (c) Isosurface of localized vacancy level illus¬ 
trating that it is composed of 5d states derived from the three 
nearest cations, (d) Isosurface of localized hole level for lowest 
energy Vx-center (e-g) Gonflgurations of the Geua —Srua —Ter 
triple complex, in which both Sr and Ge are within the first 
neighbor shell of the Br vacancy, (e) in-plane Ceua — Vbt and 
out-of-plane SrLa — Vbt- (f) out-of-plane Geua — Per and out- 
of-plane Srca — Vbt- (g) out-of-plane Geua — Vbt and in-plane 
SrLa — Vbt- 


Using the pSIC method,^® we have identified several 
self-trapped hole configurations (also known as Vk cen¬ 
ters). Analogously to the classic Vk centers in Nal, the 
self-trapped polaron configurations in LaBra involve a 
dimerization of two halide ions. We find that the Br~ 
Br distance of the most energetically favorable Vk cen¬ 
ter is 2.98 A compared to 3.6 A in the ideal lattice, see 
Fig. 4(d), with a binding/formation energy of —0.3eV. A 
more detailed exposition on polaron binding energies and 
migration barriers in LaBra will be published elsewhere. 

Knowledge of the formation energies in combination 
with the charge neutrality condition, see Sect. II, allows 
one to compute defect concentrations as a function of 
the chemical boundary conditions as done previously in 
Ref. 14. Here, it is exemplarily illustrated in Fig. 5(a), 
which shows the dependence of defect concentrations and 
free charge densities on the relative chemical potential 
of Br, A/iBr, assuming full equilibrium at a tempera¬ 
ture of 600 K. In this representation, A^Br —>■ 0 cor¬ 
responds to Br-rich conditions, implying that the ma¬ 
terial is equilibrated with respect to a Br-rich reservoir 
such as Br 2 gas, see Appendix A. In the La-rich limit, 
A/iBr —^ |-ff/(LaBr 3 ), which is equivalent to A/iLa —>■ 0. 

According to Fig. 5(a) for the widest range of chemi¬ 
cal conditions charge equilibrium is accomplished by La 
and Br vacancies, which act as acceptors and donors, re¬ 
spectively, or in Kroger-Vink notation \V{^2[ ~ ^[UBr]- 
the extreme La-rich limit the present calculations sug¬ 
gest that the La vacancies in this balance are replaced 
by free electrons, albeit at a small concentration. As the 
band gap is underestimated with respect to experiment 
(5.3 vs 5.9eV, compare Sect. II), the concentrations of 
free charge carriers are somewhat overestimated relative 
to defects. Being aware of this shortcoming we focus on 
the Br-rich limit from here on. The general conclusions 
drawn from our results are, however, entirely unaffected 
by this issue. 


B. Strontium 

After the basic properties of intrinsic defects have been 
established, one can now explore the effect of Sr incor¬ 
poration. The energetics of interstitial as well as substi¬ 
tutional defect configurations were considered as shown 
in Fig. I(c,d), from which SrLa emerges as the domi¬ 
nant form. This defect acts as a singly charged acceptor 
Sr^a over the widest range of electron chemical poten¬ 
tials with an equilibrium transition level less than 0.2 eV 
above the valence band maximum (VBM). Both the posi¬ 
tioning of the transition level in the vicinity of the band 
edge and vanishingly small structural changes between 
charge states q = —1 and 0 indicate that the defect is 
electronically shallow. 

Under certain chemical conditions interstitial Br, 
which acts as a shallow donor, can also assume low for¬ 
mation energies. From a more detailed analysis of de¬ 
fect concentrations this defect is, however, found to occur 
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FIG. 5. Defect and free charge carrier concentrations in (a) pure and (b) Sr-doped LaBrs as a function of chemical potential 
assuming full equilibrium at 600 K (also compare Ref. 14). (c) Concentrations of Br (solid lines) and La (dashed lines) vacancies 
as a function of temperature for different levels of Sr doping and assuming defect freezing at BOOK. Note that incorporation of 
Sr leads to an increase in the Br vacancy concentration by several orders of magnitude. 


generally in much smaller concentrations than Srua- 

The effect of Sr doping on intrinsic defect concentra¬ 
tions is exemplarily shown for a temperature of 600 K 
in Fig. 5(b). The vast majority of Sr is incorporated as 
Srpa and as an acceptor is balanced by Br vacancies, i.e. 
[Sr] « [SrL,j] « [V^,.]. As typical doping concentrations 
of Sr are between 50 and 200 ppm^^ this implies that 
co-doping with Sr leads to a substantial increase in the 
concentration of Br vacancy compared to pure LaBra. 

Due to the limited mobility of atoms (and thus de¬ 
fects) at low temperatures crystalline materials are usu¬ 
ally not in full defect equilibrium at low temperatures, 
say near room temperature. Rather defect concentra¬ 
tions are “frozen in” as the material is cooled down after 
manufacturing. A representative “freezing” temperature 
of 600K was assumed to generate Fig. 5(c), which shows 
the evolution of the concentrations of intrinsic defects 
both for pure and doped materials. The data demon¬ 
strates that Sr doping can be expected to increase the Br 
vacancy concentration by up to five orders of magnitude 
compared to (nominally) pure material. 

The opposite charge states of Vg,. and SrLg^ cause mu¬ 
tual attraction as already shown in Ref. 14. The interac¬ 
tion strength is quantified in Fig. 6(a) revealing a binding 
energy of —0.3eV for the nearest neighbor (Sruaber)^ 
complex. A closer inspection of the electronic structure 
of the complex reveals that the defect levels associated 
with the Br vacancy are shifted closer to the CBM by up 
to 0.4 eV compared to the isolated vacancy, see Fig. 6(b). 
The shift of the defect level can be rationalized by con¬ 
sidering that SrLg^ (unlike for example CeLj^) introduces 
a point charge-like electrostatic potential that shifts the 
local energy scale, see Fig. 6(c). Localized states such as 
the Vbi defect level are sensitive to this shift whereas the 
delocalized states that make up the valence and conduc¬ 
tion bands are unaffected, causing an effective upward 
shift of the vacancy level. 


It should be noted that there are two distinct first near¬ 
est neighbor Srua — ker configurations, which correspond 
to the Sr vacancy being oriented out-of-plane (point A) 
and in-plane (point B) relative to each other, respectively, 
compare Fig. 4. In spite of the very similar separation 
between vacancy and Sr, the out-of-plane complex has 
an equilibrium transition level that is 0.2 eV closer to 
the CBM than in the in-plane geometry. This difference 
arises from the local constraints on relaxation that affect 
the two configurations differently. A similar effect is also 
observed in connection with the Stokes shifts associated 
with Ceua-vacancy complexes, see Sect. HID. 

In short, the incorporation of Sr in LaBr 3 thus (i) in¬ 
creases the Br vacancy concentration by several orders 
of magnitude and (ii) reduces the separation between va¬ 
cancy level and CBM. These two effects have important 
implications with respect to understanding the improved 
scintillation response of Ce/Sr co-doped LaBra as will be 
discussed in detail in Sect. IV. Another important factor 
in this regard that should be mentioned here is the very 
small lattice distortion that occurs upon Sr incorpora¬ 
tion, which will be revisited in Sect. HIE. 


C. Cerium 

We now move on to consider the incorporation of Ce. 
Substitutional Ceua is energetically clearly preferred over 
Ce interstitial and Cesr configurations. It remains neu¬ 
tral over the widest range of /Ze corresponding to Ce^’*’- 
4/^5d°. If fie drops below ctr = 0.37eV the 4/ level 
is depopulated leading to a configuration that corre¬ 
sponds to Ce^~^-4f^5dP. This range of electron chemi¬ 
cal potentials is, however, typically not observed in ex¬ 
perimental settings as demonstrated by the absence of 
Ce^^ signatures.Under La-rich conditions the forma¬ 
tion energy difference between interstitial and substitu- 
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Distance to defect atom (A) 

FIG. 6. (a) Binding energies of Sr with various defects as a 
function of separation. The Srua — Vbi — Ceca triple complex 
is composed of a nearest neighbor in-plane SrLa — Vbt defect 
and CeLa at various distances, (b) Position of equilibrium 
transition and trapping levels of Br vacancy as a function 
of separation from Srua (also compare Ref. 14). The arrows 
on the right-hand side indicate the levels obtained for the 
free Br vacancy, compare Fig. 2. The data points marked 
(A) and (B) correspond to first nearest neighbor out-of-plane 
and in-plane configurations, respectively, compare Fig. 4. (c) 
Potential induced by Ce^g^ and SrL,^ calculated by spherically 
averaging the difference in local potential between defect and 
ideal cell. The orange line is a fit to 1/r. 


tional Ce is reduced but CeLa remains the most domi¬ 
nant defect. The formation energy for Ce^a is —0.71eV 
assuming equilibrium with CeBra regardless of chemical 
condition. The negative formation energy indicates full 
solubility of Ce in LaBra, which is compatible with the 
large amounts of Ce that are routinely substituted into 
the material. 

As shown in Fig. 6(a), the interaction of Ceua with 



FIG. 7. (a) Configuration coordinate diagram illustrating 

the relation between ground and excited state of Ceua in 
its neutral charge state. The ground and excited state PES 
correspond to the electronic conhgurations Ce^^-4/^5d° and 
Ce^^-4/°5d^, respectively. The spin-orbit splitting of the 
ground state PES is represented by filled and open circles, 
(b) Schematic representation of the displacements associated 
with the conhgurational coordinate. The primary change per¬ 
tains to the six out-of-plane Br neighbors of the Ce site, which 
relax inward by about 0.1 A going from the ground to the 
excited state ionic configuration while the in-plane neighbor 
positions are almost unaffected. 


other defects is very weak. Since CeLa is furthermore neu¬ 
tral under all relevant conditions, its effect on the charge 
neutrality condition and the concentrations of other de¬ 
fects is negligible. 

Upon high energy radiation or absorption above the 
band gap, Ce can be excited corresponding to the tran¬ 
sition from Ce^+-4/i5d° to Ce^+-4/°5db The latter is 
associated with the emergence of electronic levels close 
to the CBM. They are predominantly of Ce-5d charac¬ 
ter and strongly hybridized with the neighboring La-5(i 
states. The calculated ground and excited state PES for 
CeLa are shown in the configuration coordinate diagram 
Fig. 7(a). While in the ground state configuration the 
ionic positions are almost unchanged compared to the 
perfect lattice, in the excited state the nearest neighbor 
Br atoms move inward by about 0.1,A. 

We obtain excitation and emission energies of 3.56 and 
3.13/2.78 eV, where the latter two values correspond to 
final states '^F ^/2 and ‘^Ft/ 2 , respectively, see Table I. 
These data underestimate the experimental values of 
4.03 and 3.47/3.19 eV (site I in Table I) by 0.47 and 
0.34/0.41 eV, respectively, which is expected given the 
well known band gap error of DFT. The error is, however, 
systematic and affects all transitions considered here in 
approximately the same way. The difference between the 
excitation and the larger emission energy gives the Stokes 
shift, for which the calculations yield 0.43 eV in good 
agreement with the experimental value of 0.55 eV. The 
relation between the ground and excited state landscapes 
is further illustrated in the configuration coordinate dia¬ 
gram Fig. 7. 

The excitation and emission lines as well as the Stokes 

























TABLE I. Comparison of calculated and experimental data for Ce absorption and emission (4/^5d° -O- 4f°5d}). Two values are 
given in the emission column corresponding to final states of ^-FV /2 and respectively. Experimental data from Ref. 13. 

IP: in-plane relative to Vbi', OP: out-of-plane relative to Vbt- AE: total energy difference; AE'^f^: relaxation energy on excited 
state PES; AE^i^-. relaxation energy on ground state PES; arc: ionic relaxation along the reaction path. 
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CeLa 



3.56 

2.78 / 3.13 

0.43 

0.19 

0.24 

0.35 

SrLa — Vbt — CeLa 









CeLa IP, SrLa OP 

4(e) 

0.00 

3.24 

2.73 / 3.08 

0.16 

0.07 

0.09 

0.24 

CeLa OP, SrLa OP 

4(f) 

0.02 

3.33 

2.73 / 3.08 

0.25 

0.11 

0.14 

0.34 

CeLa OP, SrLa IP 

4(g) 

0.09 

3.28 

2.55 / 2.90 

0.38 

0.15 

0.23 

0.49 

Experiment (Ref. 13) 









I 



4.03 

3.47 / 3.19 

0.56 




II 



3.59 

3.36 / 3.10 

0.24 




III 



3.47 

3.27 / 3.10 

0.21 





shift for Ce in LaBra were also calculated by Andriessen 
and co-workers,^"* who obtained a value of 0.42 eV for the 
latter in close agreement with the present calculations in 
spite of differences in the computational approach regard¬ 
ing both pseudopotentials, the treatment of 4/-states, 
and the description of the excited state. In contrast to 
those calculations, however, we do not obtain an asym¬ 
metric relaxation pattern for the excited state, in which 
the Ce ion is displaced along one of the Ce-Br bonds, 
described in Ref. 34 and interpreted as a pseudo Jahn- 
Teller distortion. Rather we obtain a symmetric displace¬ 
ment pattern as described above even if the structural 
optimization is started from the asymmetric structure 
from Ref. 34. We conjecture that this discrepancy is re¬ 
lated to the lack of DFT-l-17 correction terms in Ref. 34, 
which lead to both an erroneous ordering of La-4/ and 
5d states^* and partial occupancy of the Ce-4/ level in 
the Ce^'''-4/*5d° configuration. 


D. Optical signatures of Ce complexes 

Introduction of Sr in the lattice is associated with the 
emergence of additional features in the optical spectra,*^ 
indicating the existence of at least two additional Ce sites 
characterized by different Stokes shifts and also differ¬ 
ent absorption and emission wavelengths. As discussed 
above, Ce does not exhibit a propensity to form stable 
defect clusters. On the other hand, we have shown that 
SrLa and Ver have a strong tendency to bind in two sep¬ 
arate configurations. The Ce dopant level in detector 
material can be as high as 5%. Even without defect- 
defect interactions, this implies that statistically there is 
a large probability of 14% for a Ce to be in the immedi¬ 
ate vicinity of a Srua — Ter cluster. By inspection we find 
three different nearest neighbor Ceua — SrLa — Vbt triple 
clusters, depicted in Fig. 4 (e-g). Each of these is char¬ 
acterized by the position [out-of-plane (OP) or in-plane 
(IP)] of the Br vacancy relative to Sr and Ce. 


For the Ceua IP, SrLa OP cluster, see Fig. 4(e), 
we obtain excitation and emission energies of 3.24 and 
3.08/2.73 eV. This corresponds to a Stokes shift of 
0.16eV, much smaller than the predicted value for CeLa 
(see Table I). Similarly, the Stokes shifts of CeLa OP, 
SrLa OP [Fig. 4(d)] and CeLa OP, SrLa IP [Fig- 4(e)] are 
0.25 and 0.38 eV, respectively. In the nomenclature of 
Ref. 13 we thus tentatively assign the hrst cluster to site 
III, and the remaining two to site II. 

The smaller Stokes shifts for the triple clusters can be 
rationalized in terms of smaller structural relaxations in 
the excited 4/°5d* state as measured by apc of Eq. (2), 
see Table I. We indeed hnd a direct correlation between 
the size of the Stokes shift and the amount of relaxation. 
In effect, this confirms the hypothesis of Refs. 12 that the 
Ce excitation has a weaker influence on the geometrical 
structure if another defect is in its vicinity. 


E. Solubility analysis 

The detailed investigation of Sr related defects in 
Sects. IIIB and HID was motivated by the improvement 
of energy resolution observed for LaBrsiCe co-doped with 
Sr. Similar effects were observed with Ca and Ba whereas 
doping with Li, Na, or Mg does not improve the scintilla¬ 
tion response and can even be detrimental.^^ To resolve 
these observations, the investigation of dopant related de¬ 
fects was extended to cover both the alkaline (Li, Na, K, 
Rb, Cs) and alkaline earth (Be, Mg, Ca, Sr, Ba) groups. 
It comprised the same configurations and charge states 
that were already described in Sect. IIIB. Defect for¬ 
mation energies were computed with respect to the re¬ 
spective bromide compounds (see Appendix B, Br-rich 
conditions), which are commonly used to introduce the 
dopants in the synthesis. 

The results with regard to the dependence of the 
electron chemical potential are qualitatively similar to 
Fig. l(c,d) for all dopants considered here. The forma- 
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FIG. 8. Solubility of selected dopants computed using 
Eq. (Al) and the formation energies for substitution on La. 
Li and Cs demark the upper and lower limit of solubilities for 
alkaline doping. 


tion energies of different dopants shift, however, relative 
to each other. Under Br-rich conditions the Fermi level 
is located in the lower part of the band gap. Under 
these circumstances, substitution on La is energetically 
the most favorable form for Ca, Sr, and Ba, whereas in¬ 
terstitial configurations are preferred for the Be and Mg. 
For the alkaline metals, substitution on La sites is pre¬ 
ferred but the formation energies are noticeably larger 
than for the alkaline earth metals. The solubility can be 
obtained from the formation energies for substitution on 
La and Eq. (Al), which yields the data shown in Fig. 8. 
It is apparent that the solubilities of alkaline metals are 
generally several orders of magnitude smaller than for the 
alkaline earth metals. 

Experimentally, only Ba, Sr, and Ca are found to 
achieve an improvement in scintillator response. This 
behavior can be understood on the basis of these results 
as follows: While the alkaline metals appear to be elec¬ 
tronically and geometrically well matched to the LaBra 
lattice (small distortions, shallow levels), their solubili¬ 
ties are too small to accomplish thermodynamically sig¬ 
nificant incorporation. Be and Mg exhibit small forma¬ 
tion energies and thus relatively higher solubilities. They 
are, however, associated with substantial lattice distor¬ 
tions and interstitial-like configurations, which gives rise 
to carrier scattering and an overall degradation of elec¬ 
tronic conduction. In the end, only the heavy alkaline 
earth metals (Ca, Sr, Ba) combine excellent electronic 
and geometric match (small strains, shallow levels) with 
sufficiently large solubilities, providing a rationale for the 
success of these dopants. 


IV. DISCUSSION AND CONCLUSIONS 

To summarize, a comprehensive investigation of intrin¬ 
sic and extrinsic defects in LaBrs has been conducted 


on the basis of first principles calculations within den¬ 
sity functional theory. The vacancies, Vbv and Upa, were 
found to be the most dominant intrinsic donor and accep¬ 
tor defects, respectively. The equilibrium concentration 
of Uer in the nominally pure (undoped) case is about 
10 ^^cm“^, assuming a freezing-in temperature of 600K. 
The Br vacancy is associated with deep equilibrium tran¬ 
sition, trapping, and detrapping levels located several 
tenths of an eV below the CBM. 

Sr substitutes for La and acts as a shallow acceptor. 
Assuming a Sr dopant concentration of 200 ppm, charge 
compensation to ensure overall neutrality increases the 
Vbt concentration by up to five orders of magnitude. The 
coulombic attraction between these two species causes 
formation of stable SrLa — lAr complexes with a bind¬ 
ing energy of about —O.SeV. Upon this reaction the de¬ 
fect level associated with the vacancy is shifted by up to 
0.4 eV toward the CBM. Thus, incorporation of Sr intro¬ 
duces a multitude of shallow electron traps. 

Cerium substitutes for La in a local 3-1- charge state 
and the calculations support the experimental conclusion 
that no appreciable amount of Ce 4+ is present.Since 
Ce substitutes another triply ionized atom it has no elec¬ 
trostatic interactions with other defects and hence does 
not prefer to bind to either SrLa, Vbt, or their complexes. 
The high Ce concentrations used in practice, however, 
imply that there is very high probability of 14% for a Ce 
atom to be in the first neighbor shell of a SrLa — %Br com¬ 
plex. This picture is supported by the good agreement 
between experimentally observed and computed Stokes 
shifts for the 4/ o 5d transition of isolated as well as 
complexed Ce. 

In the Ce 4f^5dP configuration, the empty 5d levels 
reside inside the conduction band while the occupied 4/ 
state is associated with a deep level inside the gap. Con¬ 
versely, in the excited 4f°5d^ configuration, an occupied 
5d level is present inside the band gap. Excluding the 
possibility of energy transfer from an exciton, we note 
that the Ce is most likely to be excited via a sequen¬ 
tial hole and electron capture. In the Ce-only doped 
case, which is almost identical to the nominally pure ma¬ 
terial, no substantial amount of electron traps can be 
expected. Eurthermore, since the Ce 4/ level is very 
deep, it is natural to assume that the initial hole cap¬ 
ture is the rate limiting step and thus that there is no 
fast mechanism to reduce the electron/hole density. This 
certainly favors Auger recombination, the rate of which 
has a cubic dependence on the excitation density and 
which has been shown to be a major quenching chan¬ 
nel at the relevant time scales for halide scintillators in 
0-scan experiments. 

Conversely, by co-doping with Sr, the electron trap 
levels not only become more shallow, which allows for 
faster trapping/detrapping rates, but the trap density 
increases by several orders of magnitude in the form of 
SrLa — Vbt complexes. If these traps are active during the 
initial thermalization stage (2-10 ps in halide systems^^) 
they will effectively reduce the free electron density. As 
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a result, a larger number of holes will remain available 
for ionization of cerium activators as the probability for 
quenching of electron-hole pairs via the Auger mechanism 
should be greatly reduced. Recent picosecond optical ab¬ 
sorption experiments have shown that energy transfer to 
europium activators in Srl2:Eu may be as fast as 400 fs,^® 
which demonstrates that very fast capture is indeed pos¬ 
sible. As each defect complex will be in close proximity 
to a Ce atom, once any of the nearby Ce atoms cap¬ 
tures a hole, Coulombic attraction serves as a driving 
force for transferring the electron from the complex to 
the activator. This suggests that non-linear quenching is 
reduced at the cost of longer decay-times. In fact, two of 
the three cerium sites discussed in Ref. 12 are associated 
with very long decay times ranging from 60 to 2500 ns 
while accounting for 20-45% of the total light output.^® 

The mechanism outlined above demonstrates that co¬ 
doping of wide gap materials, in particular scintillators, 
provides an efficient means for managing charge carrier 
populations under out-of-equilibrium conditions. In the 
case of LaBr3:Ce,Sr the co-dopant manipulates not only 
the concentrations but also the electronic properties of 
intrinsic defects (specifically Ver) without introducing 
additional gap levels. This leads to the presence of shal¬ 
low electron traps that can localize charge carriers on 
a nanosecond time scale, effectively deactivating charge 
carrier recombination channels. The principles of this 
mechanism are therefore not specific to the material con¬ 
sidered here but can be adapted for controlling charge 
carrier recombination in other wide gap materials. 
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Appendix A: Defect thermodynamics 

In the dilute limit the equilibrium concentration of a 
defect depends on its free energy of formation, AG/, via 

c = coexp (-ACz/fc^T), (Al) 

where Cq is the concentration of possible defect sites. The 
formation free energy AG/ is usually approximated by 
the formation energy AE/, which is legitimate if the 
vibrational entropy and the pressure-volume term are 
small.®® The formation energy AE/ of a defect in charge 
state q is given 

AEf = Edet — Aid + ^(eVBM + Me) — ^ (A2) 

i 

where Adef is the total energy of the system containing 
the defect and Aid is the total energy of the ideal host. 
The second term describes the dependence on the elec¬ 
tron chemical potential, Me, which is measured with re¬ 
spect to the valence band maximum (VBM), eyBM- The 
formation energy also depends on the chemical potentials 
of the constituents as given by the last term, where the 
difference between the number of atoms of type i in the 
ideal cell with respect to the defect cell is denoted by 
Ani- The chemical potential Mi of constituent i can be 
rewritten as Mi = Mi“*^ + AMi where Mi“^* denotes the 
chemical potential of the standard reference state. Ne¬ 
glecting entropic contributions the chemical potentials of 
the reference phases can be replaced by their cohesive 
energies at zero Kelvin. Note that the summation in 
the last term also includes dopant or impurity species, 
whence one has to consider the source of the dopant or 
impurity atom when discussing formation energies and 
solubilities (see Sect. HIE). The chemical potentials of 
La and Br are coupled to each other via the formation 
enthalpy of the compound, i.e., A/iua + 3A/tBr = AiJ/. 
This implies that specifying either A/iLa or A/iBr is suffi¬ 
cient to determine the chemical equilibrium with respect 
to the host. Following common practice we refer to La 
and Br-rich conditions, which correspond to A/tua = OeV 
and A/tBr = OeV, respectively. 

Appendix B: Thermodynamic boundary conditions 

When calculating defect formation energies of intrinsic 
defects according to Eq. (A2) it is sufficient to consider 
the chemical potentials for La and Br only as indicated 
above; specifically, for La-rich conditions 

AMLa = 0 and A/iBr = A/(LaBr3)/3 (Bl) 

whereas for Br-rich conditions 

A/iBr = 0 and A/iLa = A/(LaBr3). (B2) 

Once extrinsic elements have to be taken into account 
additional conditions must be invoked. To be specihc 
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consider the case of Sr, which introduces one additional 
chemical potential, fisr = Msr + ^MSr, in Eq. (A2). One 
could assume the dopant to be in equilibrium with its 
elemental (metallic) form, which implies 

Ansr = 0. (B3) 

Sr and Br can, however, react to form SrBr2 and it is 
therefore more appropriate to consider the equilibrium 
with respect to the compound, which is also used exper¬ 
imentally for introducing the dopant during synthesis, 

A^iSr + 2A/rBr = ///(SrBra). (B4) 

Since this reaction involves an element of the host the 
intrinsic boundary conditions (La and Br-rich) explicitly 
affect the condition for the chemical potential of Sr. In 
Br-rich conditions 

A^Br = 0 

AfiSr = iL/(SrBr2) (B5) 

while in the La-rich limit 

AfiBr = ^///(LaBra) 

^AAiSr = Lf/(SrBr2)-^i//(LaBr3). (B6) 

The extension to other elements is straightforward. For 
the chemical potentials of the alkaline metals for example 
one obtains 

A^Na = iL/(NaBr) Br-rich (B7) 

A^Na = iL/(NaBr) — ^iJ/(LaBr3) La-rich. (B8) 

These conditions are used in Sect. HIE to determine the 
solubilities of various dopants in LaBra. 


Appendix C: Finite-size scaling 

Given the small dielectric constant of LaBra and the 
large defect charge states that need to be considered 
it is crucial to properly correct for both electrostatic 
and elastic image charge interactions. Various correc¬ 
tion schemes have been proposed for this purpose but 
ambiguities remain.In the present work we therefore 
resort to finite-size scaling, which if computationally af¬ 
fordable is expected to yield the most reliable results. 
Finite-size scaling is most commonly based on “simple” 
multiples of the underlying primitive cell. For example 
for diamond and zincblende structures supercells based 
on simple cubic, body-centered cubic, and face-centered 
cubic unit cells are often used.^®’^®’"*^ Since each of these 
cells is associated with a different Madelung constant the 
size dependence of for example the monopole-monopole 
correction, which is the leading electrostatic interaction 


term,"*^ will differ between these cells. It is therefore ad¬ 
vantageous to consider scaling among a set of self-similar 
cells. 

The direct application of this principle to the hexag¬ 
onal unit cell of LaBra would allow the construction of 
only a very small number of supercells, which in turn 
would limit the reliability of the finite-size scaling proce¬ 
dure. To resolve this situation we devised a systematic 
way to construct “optimal” supercells. Optimality here 
implies that we seek to find supercells that for a given 
size (number of atoms) most closely approximate a cu¬ 
bic cell shape. This approach ensures that the defect 
separation is large and that the electrostatic interactions 
exhibit a systematic scaling. (Recall that for example the 
monopole-monopole interaction is given by the Madelung 
constant, which is only dependent on the shape of the 
unit cell). 

The cubic cell metric for a given volume O is 

hcub = (Cl) 

which in general does not satisfy the crystallographic 

boundary conditions. The l2-norm provides a convenient 
measure of the deviation of any other cell metric from a 
cubic shape (“acubicity”) 

A,(h) = ||h-heub||2 . (C2) 

Cell metrics that are compatible with the crystal symme¬ 
try can be written as integer multiples of the underlying 
primitive unit cell hp, i.e. 

h = Php where P (C3) 

The optimal cell shape multiplier for a given cell size is 
then obtained as follows 

Popt = argnnn{Ac(Php)| detP = Nuc} , (C4) 

where N^c is the desired system size in multiples of the 
primitive unit. This approach is generally applicable and 
can be readily generalized to optimize toward other su¬ 
percell shapes, e.g., face-centered or body-centered cubic. 

A series of supercells was generated based on Eq. (C4) 
for all possible sizes up to 752 atoms {Nuc = 94). The 
lowest values of Ac achieved in this fashion are shown in 
Fig. 9. We empirically find that supercells, for which Ac 
is lower than 0.02, are sufficiently close to a cubic shape 
for our purposes. The sizes for which this limit is reached 
are indicated by the filled red symbols in Fig. 9. 

Using supercells containing between 24 and 496 atoms 
defect calculations were carried out for Br and vacan¬ 
cies and antisites in charge states that were identified as 
relevant ones on the basis of earlier 96-atom cell calcu¬ 
lations. The thus obtained configurations were analyzed 
as follows. Formation energies were computed using the 
thermodynamic formalism described in Appendix A. The 
formation energies obtained according to Eq. (A2) are re¬ 
ferred as “raw” data in the following. In addition we con¬ 
sidered the effect of potential alignment (PA) and image 
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FIG. 9. Lowest deviations from a cubic shape obtained for 
different cell sizes via Eq. (C4) and Eq. (C2). 

charge corrections. The potential alignment correction 
amounts to a term qvpA- In the present work we de¬ 
termined wpA by taking the difference of the electrostatic 
potential between defect and ideal supercell, with the po¬ 
tential being measured by a test charge at the ionic site 
farthest from the defect center. For the image charge cor¬ 
rection we adopted the simplified correction described in 
Ref. 25, which involves the addition of a term q'^'^\Ejnp\ 
where E^p is the electrostatic energy associated with a 
periodic array of point charges according to the super¬ 
cell metric taking into account dielectric screening due to 
both electrons and ions. This correction term should re¬ 
flect both monop ole-monop ole and monopole-quadrupole 
terms. The resulting expression for the formation energy 
is 

^Ef = AF/j -|- qvpA + 9^ 21-E'mp I • (Cb) 

The leading terms in Eq. (C5) should scale with 
and V (or and N), therefore we also fit formation 

energies obtained from Eq. (A2) to the following expres¬ 
sion 

AEf = AEf{N) + a A" -b 6A, (C6) 

where AE^, a and b were treated as fit parameters. 

The results of this analysis are exemplified in Eig. 10, 
which illustrates the scaling for La vacancies in various 
charge states. Analysis of e.g., Br vacancies, antisites. 


and interstitials yield very similar plots. In all cases we 
find that the combination of potential alignment and im¬ 
age charge corrections yields formation energies that ap¬ 
proximate the infinite limit rather well. The accuracy of 
the thus corrected data is, however, limited by the accu¬ 
racy associated with the determination of the potential 
alignment correction and the approximitive nature of the 
image charge correction. Note that the corrections do not 
account for elastic image interactions. The comparison 
of extrapolated values, which do include elastic effects, 
and corrected data, which do not, indicates that elastic 
interactions are, however, negligible in the present case. 
This conclusion is further supported by the fact that the 
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FIG. 10. Finite-size scaling of formation energies of La vacan¬ 
cies in different charge states. Open symbols refer to forma¬ 
tion energies calculated via Eq. (A2) whereas filled symbols 
belong to formation energies subjected to potential alignment 
and image charge corrections according to Eq. (C5). Solid 
lines are fits to Eq. (C6) and dashed lines indicate the av¬ 
erage over the corrected formation energies. The numbers 
next to the lines indicate the formation energies obtained by 
extrapolation of the uncorrected data. In the case of the neu¬ 
tral vacancy (q = 0) there are no potential alignment and 
image charge corrections. Formation energies have been com¬ 
puted for Br-rich conditions and an electron chemical poten¬ 
tial pe = 0.6 eV, which was chosen to obtain a visual separa¬ 
tion of the different charge states. 


formation energies of neutral defects are almost indepen¬ 
dent of system size for N > 96. 
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